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ABSTRACT 

We compute rates of period change (P's) for the 215s mode in G117-B15A and 
the 213s mode in R548, first for models without axions, and then for models with 
axions of increasing mass. We use the asteroseismological models for G117-B15A 
and R548 we derived in an earlier publication . For G117-B15A, we consider two 
families of solutions, one with relatively thick hydrogen layers and one with thin 
hydrogen layers. Given the region of parameter space occupied by our models, 
we estimate error bars on the calculated P's using Monte Carlo simulations. 
Together with the observed P for G117-B15A, our analysis yields strong limits 
on the DFSZ axion mass. Our thin hydrogen solutions place an upper limit of 
13.5 meV on the axion, while our thick hydrogen solutions relaxes that limit to 
26.5 meV. 

Subject headings: Dense Matter — Stars: Oscillations — Stars: Variables: Other 
— Stars: White Dwarfs — Elementary Particles 



1. Astrophysical context 

G117-B15A and R548 are pulsating white dwarfs with atmospheres dominated by hy- 
drogen. These stars are called DAVs or ZZ Ceti stars. Their effective temperatures are 
close to 12,000 K. They are non-radial, g-mode pulsators (the restoring force is buoyancy). 
G117-B15A is the most stable optical clock known. The rate of change in pe riod (P= dP/dt) 



of its most stable mode (215s) is 3.57±0.82 x 10"^^ s/s flKepler et al.ll2005h . That is, G117- 
B15A loses one "tick" every 1.7 billion years. We can use this well-measured P to constrain 
among other things the mass of the axion, the best motivated dark matter candidate and in 
general, the emission rate of weakly interacting particles. While measuring a P for R548's 
most stable mode (near 213s) is less trivial because it is a doublet, the current upper limit. 
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5.5 ±1.9 X 10 s/s (iMukadam et al.ll2003l ). also indicates that in the near future R548 may 
become useful to constrain the axion mass as well. 

A P of order 10~^^ for DAVs is consistent with a slow increase in period due to the 
cooling of the interior. It is an evolutionary timescale — a measure of how fast the star is 
cooling. Knowing P, we can determine how much energy a star like G117-B15A is losing 
every second . This energy loss includes of course the radiation of photons (Mestel cooling, 
Mestellll952l ). but in principle could also include the emission of weakly interacting particles. 
While neutrinos immediately come to mind, the Standard Model of Part icle Physics predicts 



that neutr ino emission should be negligible in G117-B15A and R548 (jWinget et al.l 12004 



Kinj 120071 ). One possibility includes axions and we apply our method to the latter in the 
present study. 



Turner 


1990; 


Kim 


2007) 



particles, axions are currently favored candidates for dark matter. But they have not been 
discovered (neither have supersymmetric particles) and the theory of axions fails to place 
any constraint on their mass. The possible contribution of axions to dark matter depends 
on their mass. The mass of axions determines how strongly they interact with the matter 
we know, with more massive axions interacting more strongly. In turn, the stronger the 
interaction of axions with matter or light, the larger their emission rate. With pulsating 
white dwarfs, we can constrain the axion emission rates and therefore their mass. 

We determine limits on the emission of weakly interacting particles by comparing their 
cooling effect on a white dwarf model with the observed cooling rates of that star (given by 
P). The strength of those limits depends on the uncertainties involv ed in the rn e asure ment 
of the P and i n the modeling. This method has been used before by llsern et al.l (119921 ) and 
Corsico et al.l (120011 ) to derive an upper mass limit for axions. While the limits on the axion 
obtained b y those authors were stronger than those obtained through other methods (see 
KimI 120071 for a review across fields of study), we now have the tools to greatly strengthen 
them. We present here the first study that takes into account the modeling uncertainties in 
a systematic way, allowing us to derive a strong upper mass limit for the axion. We present 
our results in a way that facilitates their application to other weakly interacting particles. 
We also publish for the first time the P we would expect to observe in R548 as a function of 
the emission rate of weakly interacting particles. Together with a future measurement of P 
for R548, these results can be used to reinforce the limits given by G117-B15A. 



In iBischoff-Kim et al.l (120071 ). thereafter paper I, we used a systematic, fine grid search 
to find best fit models to G117-B15A and R548. We were able to define the regions of 
parameter space the fits occupied. In this work, we use Monte Carlo simulations to calculate 
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P's that include theoretical uncertainties, for varying axion mass. We begin with a brief 
introduction to a xions in section 2. I n sec tion 3, we give an overview of the work done by 



Isern et al.l (Il992l ) and lCorsico et al.l (120011 ) . We detail our method in section 4, discuss our 



results in section 5, and conclude in section 6. 



An introduction to axions 



The Standard Model of particle physics does not explain why Charge-Parity (CP) should 
be violated in weak interactions but not in strong interactions. One would therefore expect 
the neutron to have an electric dipole moment. Experiments have placed an upper limit on 
the neutron's electric dipole of 10~^^ e -cm (electron charge x length of dipole), 10 orders of 
magnitude below the predicted value (lTurnerlll990l ). 



In more technical terms, in QCD, the neutron elect ric dipole mom ent arises from a CP 
(as well as C and P) violating term in the lagrangian (jTurneiill990l ). This electric dipole 
moment is predicted to be of order 5 x 10~^^6' ecm, where ^ is a free parameter in the 
theory. The experimental limit of 10~^^ecm means that 6 is less than 10~^°. Theoretical 
physicists do not understand why 6 has to be so small. This is called the strong CP problem 
and can be solved elegantly by the introduction of a new symmetry. This symmetry, added 
to the Lagrangian of the fundamental interactions, is called the Peccei-Quinn symmetry 
(jPeccei fc Quinru 119771 ). The spontaneous breaking of this symmetry gives rise to axions. 
The theory does not place any limit on their mass. 

Axions can couple to photons, electrons (leptons) and nucleons (baryons). The relative 
strength of each coupling depends on one or more of the Peccei-Quinn charges of the u 
and d quarks and the electron (denoted Xu, Xd, and Xe respectively.) Those charges can 
be between and a number of order 1 and are not constrained by the theory. This gives 
rises to a continuura of a xion models. Two simple models include the KVSZ model ( iKim 



979 



(iDine et al 



1981 



Shifman et al.lll980f). where X^ = (no coupling to electrons), and the DFSZ model 
Zhitnitskvl[l980h . where A„ ~ A^ ~ A^ ~ 1. 



In this paper, we focus on DFSZ axions, those that interact with electrons. For a more 
co mplete review of the different axion models and their interaction with light and matter. 



see 



KimI (120071 ). In white dwarf interiors, the dominant axion emission mechanism would 



be electron bremsstrahlung, where an electron emits an axion as it gets accelerated in the 
neighborhood of an ion. The axion emission rate (per unit mass) for that process may be 
expressed as (INakagawa et al.l 1 19881 ) 



1.08 X 10^^ ergs g"^ s'^a — T^ F(T, p). 
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gl^/AiT, and gae, the strength of the axion-electron cquphn g; is 



where T7 = T/IO^K, a 
given by equation [2] below. F(T,p) is a numerical fit provided by lNakagawa et al.l (Il988l ). It 
is of order 1 throughout most of the interior of a typical white dwarf model (e.g. figured]). 



The strength of the axion-electron coupling is 

gae = (2.8 X 10 



11 maCOS^/? 



1 eV 



(2) 



where ma and f3 are free parameters. For DFSZ axions, the coupling to electrons is 10^ 
or ders of magnitude greater than the coupling to photons (assuming cos^/5 ~ 1). Following 



m 



Isern et al.l (jl992|)'s footsteps, we shall not make make any assumptions on the value of (3 
and simply state that the quantity we are constraining is ma cos^ (3. We shall, however, refer 
to it simply as the "axion mass" . 



3. Early work 



In a pioneering work, llsern et al.l (1l992l ) used G117-B15A's P t o obtain a limit on the 
axion mass. At the time, the P measured (12.0 ± 3.5 x 10~^^s/s, iKepler et al.l Il99ll ) for 
that star was uncertain, and much h igher than the one expected from simple M estel cooling. 
Using models available at t he time (WoodI 199ol : D'Antona fc Mazzitelli 1989 ) and a simple 
semi-analytical treatment, llsern et al.l found an average axion mass of 8 meV. Individual 
values, depending upon the model chosen and value of observed P considered (P — AP, P, 
and P -|- AP) allowed a range between meV and 20 meV for the axion mass. 



Corsico et al.l (1200 ll ) revisited the problem with a new mea sured yalue o f P, that had 



since decr eased to what was expected from si mple Mestel co oling ( Mestel 19521 ): (2.3 ±1.4) x 
10~^^s/s ( IKepler et al.ll2000l ). In their work. ICorsico et al.l performed an asteroseismological 
study of G117-B15A to find its mass, helium layer mass, and hydrogen layer mass. To 
reduce the number of parameters to fit, they fixed the i nternal composition to that found 
from stellar evolution calculations by ISalaris et al.l (119971 ) , and the effective te mperature to 
the latest spectroscopic estimate at the time, 11,620 K (IBergeron et al.lll995l ). Their best 
fit model had a mass of 0.55 Mq, a helium layer mass Mne = 10~^M^, and a hydrogen layer 
mass Mhe = 10~^M*. On the average, the periods of that model are 5 seconds away from 
the observed periods. In contrast, the models we use in the present study fit to better than 
1 second. 



Corsico et al.l (I2OOII ) considered only small uncertainties in effective temperature (200K) 



and found that they led to a 4% uncertainty in the calculated P's. They also considered 
larger uncertainties in mass and central oxygen abundance and found that those had less 
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of an effect on the model's P's than the effective temperature. A mass uncertainty of .02 
Mq (3%) led to a 6% uncertainty in P's. And considering a full range of core composition 
(0% carbon to pure carbon) changed the P's only by 5%. They concluded that uncertainties 
other than the one in the measured P were insignificant and could be ignored altogether. 
They found a tight constraint on the axion mass: 4 meV. If we adopt Corsico et al.'s best 
fit parameters for G117-B15A, and follow the same method using our models, we obtain a 
similar limit on the axion mass. 

To place things in prospective, the best DFSZ axion mass limits aside from that of 
Corsico et al. and the ones presented here come from the observation of Red Giants 

and Horizontal Branch stellar populations in clusters. If axions were massive enough, they 
would cool the core of stars and have an observable effect on the morphology o f those 
populations. The best limit found using this method is 9 meV (IRaffelt fc Weisslll995l ). Such 
limits depend on the assumed core mass in the models, a rather uncertain quantity. The 
pulsating white dwarf method does not suffer from such large uncertainties. 



Since the work done bv I Corsico et al.l (120011 ). a newer value of P = 3.57±0.82 x lO^^^s/s 
has been obtained (IKepler et al.ll2005l ). The new value has a smaller measurement error, 
and could help constrain the axion mass better. More importantly, we take into account 
additional modeling uncertainties that have been ignored in previous analyses. The better 
we account for those uncertainties, the stronger our limit on the axion mass will be. We will 
see that this leads us to different conclusions than the authors mentioned above. 



4. Method 

4.1. Axion emission in our models 

For a mixture of carbon and oxygen, as in the interior of our white dwarf models, the 
axion emission rate is given by 

Ca = Xcec + Xoeo, (3) 

where Xc and Xq are the carbon and oxygen abundance, and ec and epsilorio are the axion 
emission rates in pure carbon and pure oxygen respectively. Combining equation [T] and [3] we 
have: 

e, = 1.08 X 10^3 a [3X,F,{p, T) + 4X,Fo(p, T)] . (4) 

We distinguish F(T, p) in pure carbon and pure oxygen by using subscripts. 

Integrating over the mass of the model, we obtain the axion luminosity. While we carry 
the exact integration numerically in our models, it is useful at this point to write down 
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an approximate expression for the axion luminosity. We use the internal properties of a 
fiducial model, chosen because it matches G117-B15A's periods and P for the 215s mode 
very closely. This fiducial model has Tgfj = 12656 K, M^, = 0.602 Mq, a helium layer mass 
of 3.55 X 10~^M^,, and a hydrogen layer mass of 4.79 x lO^^M^.. We display the internal 
properties of the fiducial model in figure [H 




0.4 0.6 
M/M* 



0.4 0.6 

M/M* 



Fig. 1. — Internal properties of the fiducial model. Left panel: carbon and oxygen abun- 
dance. Right panel: log of the temperature in K, log of the density in g/cc, and Fc (see 
text). 
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FoUowing the example set forth in llsern et aLl (119921 ). we make the following approxi- 
mations: 1) We set T = Tccntcr = 1-2 X 10 ''K throughout the model; 2) We also set Fc and 
Fq constant (Fc ~ 1.1 and Fq ~ 1.2). The axion luminosity is then 



La = 1.08 X 10^^ aT^M, ( 3Fc / Xcdm + 4Fo/ Xodm 



(5) 



where m = Mi./M^,. For the fiducial model's core composition, dm ^ 0.4 and Xq dm % 
0.6. We obtain 

La = 1.13 X 10^^ a ergs s~^ (6) 
This approximate expression is 7% accurate for our fiducial model. 



4.2. Asteroseismological fits 

Before we include axions in our models, we need to make sure we have reliable models 
of the stars under study, G117-B15A and R548. We performed our own asteroseismological 
analysis of those stars in paper I. We varied 4 parameters, including the effective temperature 
and the helium layer mass, Mne (in addition to the stellar mass and the hydrogen layer 
mass, Mh). In paper I, we explored systematically a number of sources of uncertainties in 
our models and concluded that we could adequately account for them by treating equally 
all models that matched the observed periods to better than 1 second. With that cut- 
off, we defined the regions of parameter space our models occupied. We found that our 
best fit models were hot and/or massive compared to the spectroscopy, though they were 
not completely inconsistent with the latter. Work is in progress to resolve this apparent 
discrepancy. Our calculations are self-consistent and we proceed to derive meaningful results. 
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Fig. 2. — A sample simulated random population for G117-B15A and R548 in the M^. — Teg 
plane and log(MH) — log(MHe) plane. The small symbols show the location in parameter 
space of the simulated population, while the squares and diamonds are best fit grid models 
from paper I. 
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4.3. Monte Carlo Simulations 

Having narrowed down the regions of parameter space occupied by best fit models for 
G117-B15A and R548, the next step is to derive rates of period change for those models along 
with the associated uncertainties. This is a classic Monte Carlo problem. For G117-B15A, 
we chose T^s, AM (instead of M^), Mne, and Mh as our 4 independent variables, where AM 
is defined by 

M, = 6 - a TefT + AM, (7) 

and a and b are defined so that the simulated populations occupy the region of parameter 
space in the M* — Tefr plane defined by the results from our astereseismological fits. For 
each family of solutions we found in paper I, without including axions just yet, we generated 
a random population (N=500). We display the distribution of the models of each of the 
populations in figure [2l 

We have two possible families of fits for G117-B15A. One has a thicker hydrogen layer 
(Mh = 6.3 X 10^''M^) and occupies the high mass, low effective temperature region of the 
M,, — TefT plane. The other has a thin hydrogen layer (Mh < 6.3 x 10~'^M*) and occupies the 
higher effective temperature (lower mass) region of the M* — Tefr plane. The thin hydrogen 
fits to G117-B15A are very similar to the fits to R548. While this makes the thin hydrogen 
fits more attractive than the thick hydrogen fits (by "thick" we mean "not as thin"), we 
cannot discard the latter as a viable possibility on that basis alone. We shall consider both 
families when constraining the axion mass. 

For each model in each population, we proceeded with calculating a P. We checked that 
the period and P distribution for each of the populations were close to Gaussian and found 
that they were. For each population (family of fits) we calculated the mean period and P. 
We list the results in tabled] with one sigma error bars. Remember that these results assume 

Table 1: Average periods and P's for G117-B15A and R548, without axions. 





G117-B15A 


R548 


Mh 


6.3 X 10"^ 4 X 10"^ 


4 X 10-« 


P (observed) [s] 
P (calculated) [s] 


215.197389 
214.9 ±2.5 215.7 ±2.0 


213 
213.1 ± 1.8 


P (observed) [10"^^s/s] 
P (calculated) [IQ-^^s/s] 


3.57 ±0.82 
1.92 ±0.26 2.98 ±0.17 


< 5.5 ± 1.9 
2.91 ±0.29 
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axions do not exist and that the coohng is due entirely to photons. For G117-B15A, the 
thick hydrogen fit P is two sigma below the observed value. However, we cannot discard this 
solution based on a low P alone, as we can always increase the cooling rate and raise the 
value of P by adding axion emission to the models. 

Next we included axion production in our mode ls and repeated the procedure for each 



different axion mass. We confirm an important fact ICorsico et al.l (120011 ) first pointed out: 
the periods themselves are insensitive to the axion mass (for ma up to 30 meV at least). The 
reason for this small dependence is that axion emission does not introduce any "bumps" into 
the Brunt- Vaisala frequency. It merely slightly lowers it throughout the degenerate core. 
This allows us to meaningfully compare the calculated P's with the observed value. 



5. Results 

In figure Owe show P as a function of the axion luminosity for G117-B15A and R548. 
While we included axions specifically in our models, the P's depend only on the total lu- 
minosity, not on the exact nature of the weakly interacting particles responsible for the 
"unseen" loss of energy. Our results may therefore directly be generalized to weakly inter- 
acting particles other than axions. 

The thin hydrogen solution for G117-B15A places an upper limit of 1.1 x 10"'^^ erg/s 
for the energy loss rate due to weakly interacting particles. For axions, this translates to an 
upper limit of 3.78 x 10"^'^ on the axion-electron coupling constant gae (13.5 meV on the axion 
mass). The thick hydrogen family of fits (taken alone) constrains the "unseen" luminosity to 
be between 2.7 x 10"^^ and 4.4 x 10"^^ erg/s. These limits translate to 2.9 x 10"^^ < gae 
< 7.4 X 10^^'^ (axions between 10.4 and 26.5 meV in mass). Combining the two families 
together, we can place a conservative upper limit of 4.4 x 10^'^^ erg/s on the luminosity of 
weakly interacting particles that would be emitted in the interior of G117-B15A, and a 26.5 
meV upper limit on the axion mass. For R548, we are still awaiting a better determined P 
to allow us to obtain a potentially tighter constraint on the emission of weakly interacting 
particles in white dwarfs. 



G117-B15A 



R548 




"Unseen" Luminosity [IC^Sl erg/s] "Unseen" Luminosity [10'^31 erg/s] 



Fig. 3. — P's as a function of the luminosity due to the emission of weakly interacting 
particles (in that case, axions) for G117-B15A and R548. The observed P for G117-B15A is 
indicated by the horizontal solid line and its 1 sigma error bars by the dashed lines above 
and below it. The shaded regions inside the dashed horizontal lines indicates the range of 
calculated P's consistent with the observed value (to 1 sigma). 
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6. Conclusions 

We computed P's for the 215s mode in G117-B15A and the 213s mode in R548 first 
for models without axions, and then for models with axions of increasing mass. The thin 
hydrogen, zero axion models for G117-B15A have a P of 2.98 ± 0.17 x 10~^^s/s, consistent 
with the observed value of 3.57±0.82 x lO^^^s/s. The zero axion, thick hydrogen models have 
a P of 1.92 ± 0.26 x 10~^^s/s, 2 sigma below the observed P and require additional cooling 
through axions to become consistent with that observed value. The zero axion models for 
R548 have a P similar to the thin hydrogen P for G117-B15A (~ 3 x 10~^^s/s), consistent 
with the best current limit on P for R548's 213s mode. 

If we believe that G117-B15A must be similar to R548, then we find an upper limit 
on the axion luminosity of 1.1 x 10^'^^ erg/s, corresponding to an axion-electron coupling 
constant gae = 3.78 x 10^^'^ , or an axion mass of 13.5 meV. A second family of fits with 
thicker hydrogen layers relaxes that limit to < 4.4 x 10~^^ erg/s, gae < 7.4 x 10~^^ or 
ma < 26.5 meV. What makes this limit unique is not the fact that it is lower than previous 
limits. What makes it unique is its strength. It is based on an unbiased exploratio n of 



the parameter space of the models. Compared to typical axion searches (see iKiml 120071 for a 
review) , the present analysis involved well-known physics and relatively small, well-identified 
uncertainties, considered systematically. Improved asteroseismological models will help us 
choose between the two families of best fits for G117-B15A and improve this limit. 

An even longer data baseline for G117-B15A will reduce the size of the error bars on the 
measured P and further narrow the allowed range for the axion mass. We can soon expect 
a useful P for R548 as well, which will reinforce the results of this study. As a part of his 
planet search around pulsating whit e dwarfs, MuUa lly has determined upper limits of 10^^'^ 
s/s for an additional 9 white dw arfs (Mullallv 2007h and ha s identified more blue edge DAVs 



with stable modes. With Argos jMukadam fc Nathei][2005l ) . the high-speed photometry data 



collected at McDonald observatory as a part of the on-going observation of hot DAVs should 
begin to yield measurements of P's for those stars within the next 3 to 5 years. It is merely 
a question of time until we can place an even tighter constraint on the emission of weakly 
interacting particles in these stars. 

In the present study, we assumed that any unseen source of energy loss was due to 
axions. The results we present can readily be used to constrain the emission of other weakly 
interacting particles possibly produced inside white dwarfs, such as supersymmetric particles. 
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